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Abstract 


In this study, the pressure variation and the flow distribution in the manifold of a fuel-cell stack are simulated by a computational fluid dynamics 
(CFD) approach. Two dimensional stack model composed of 72 cells filled with porous media is constructed to evaluate pressure drop caused by 
channel flow resistance. In order to simplify this model, electrochemical reactions, heat and mass transport phenomena are ignored and air is treated 
as working fluid to investigate flow distribution in stacks. Design parameters such as the permeability of the porous media, the manifold width and 
the air feeding rate were changed to estimate uniformity of the flow distribution in the manifold. A momentum-balance theory and a pressure-drop 
model are presented to explain the physical mechanism of flow distribution. Modeling results indicate that both the channel resistance and the 
manifold width can enhance the uniformity of the flow distribution. In addition, a lower air feeding rate can also enhance the uniformity of flow 
distribution. However, excessive pressure drop is not beneficial for realistic applications of a fuel-cell stack and hence enhanced manifold width is 


a better solution for flow distribution. 
© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


A proton exchange membrane (PEM) fuel cell is an electro- 
chemical device where chemical energy is directly converted 
into electricity by using hydrogen as fuel and oxygen as oxidant 
[1]. A single PEM fuel cell can offer a voltage of ~0.6-0.7 V. In 
order to increase the voltage sufficient for practical operations, 
many single cells are serially connected to fabricate the fuel-cell 
stack. Sufficient reactant feed is necessary for each cell to attain 
high performance and stable operation of stacks. The gas mani- 
folds play a key role in the uniform distribution of reactant gases 
to the individual cells. A good manifold design is important to 
achieve the above requirement. An uneven flow distribution will 
result into performance loss as well as the stack may cease to 
function. 

Till now most of the PEM fuel cell modeling works mainly 
focus on electrochemical, heat and mass transport phenomena 
in a single cell. The single cell models reported by Bernardi 
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and Verbrugge [2,3], Springer et al. [4], Nguyen and White [5], 
Gurau et al. [6], are some of the examples. Recently, the Penn- 
sylvania State University [7—12] in the US and University of 
Victoria [13-15] in Canada are reporting impressive modeling 
job. 

The above mentioned single cell models deal with the 
construction of the fuel cell model based on the numerical tech- 
niques to simulate the related phenomena that occur within the 
fuel cell. The model parameters are normally verified with the 
experimental results. The models help to predict the perfor- 
mance characteristics and understand the physical phenomena 
that occur within the fuel cell. 

In recent years, increasing number of researchers are focus- 
ing on the modeling of fuel-cell stack. The models are mainly 
focusing on the flow distribution and water-thermal management 
of stacks. Usually, the modelers ignore the complicated issues 
like electrochemistry and transport process in order to simplify 
the models. Koh et al. [16] have reported a numerical model 
to investigate pressure variation and flow distribution of stacks. 
This model considers channels filled with porous media to eval- 
uate pressure drop caused by channel configuration. In addition, 
it compares flow distribution for different overall gas flow pat- 
terns. Karimi et al. [17] have developed a stack model with fluid 
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Nomenclature 
a link coefficient 
Fẹ wall shear force (N) 


permeability (m7) 

mass flow rate (kg s7!) 

normal direction relative to outlet ports 
number of cells 

production rate of k (m? s73) 
pressure (Pa) 

total mass flow rate of stacks (kg s7!) 
source term 

pressure action area (m?) 

time (s) 

turbulence intensity 

velocity vector (m s7!) 

velocity vector (m s7!) 

Cartesian coordinates (m) 


5 


= LAND NaO Srp, 


Greek letters 

physical value 

exchange coefficient 

tensor 

turbulent dissipation rate (m? s73) 
porosity 

turbulent kinetic energy (m? s7?) 
viscosity (kg m7! s7!) 

density (kg m~?) 

shear stress (N m7?) 
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Subscripts and superscripts 


C empirical constants 

eff effective property which accounts for porosity 
in inlets 

Lj components of a vector in Cartesian coordinates 
sys system 

t turbulent 

o empirical constants 

0) general dependent variable 


mechanics method. Reactant species such as hydrogen, oxygen, 
nitrogen, etc. were introduced into this model so that it can eval- 
uate the relationship between the consumption rate of reactants 
and stack output power. Yu et al. [18] have proposed a water- 
thermal management model to investigate PEM fuel-cell stacks 
made from Ballard power system, Canada. This model can pre- 
dict temperature, voltage, power, efficiency, etc. under steady 
and transient operations of stacks, respectively. Promislow and 
Wetton [19] proposed a steady heat transfer model of PEM fuel- 
cell stacks composed of parallel cooling channels and infinitely 
thin membrane electrode assembly (MEA). This model predicts 
not only the local temperature difference between coolant and 
MEA, but also the cell temperature variation. 

The stack modeling mentioned above is mainly concerned 
with developing a numerical model characterizing gas flow dis- 


tribution and water-thermal management. The above models 
were based on numerical analysis but not on computational fluid 
dynamics. These models could not investigate detailed stack 
flow field configuration and transport phenomena. Although 
simplified models can evaluate stack performance characteris- 
tics rapidly, it is also possible to obtain wrong results due to the 
over-simplification of the model. 

Experimental data associated on the gas or coolant flow 
distribution is usually difficult to obtain. Therefore, indirect mea- 
surements like voltage or temperature distribution are used to 
understand the physical phenomena that occur in a fuel-cell 
stack. Due to the lack of experimental data, fuel-cell stack 
numerical models only provide theoretical explanations and 
understanding about performance characteristics. In this study, 
since there is no appropriate experimental data available, phys- 
ical mechanism about flow distribution is investigated only by 
using purely theoretical approach. But it is indeed an available 
and acceptable method for studying phenomena that could not 
be experimentally measured or calculated. 

This study presents a computational fluid dynamics model of 
PEM fuel-cell stack to investigate pressure variation and flow 
distribution. A two-dimensional model consisting of 72 chan- 
nels filled with porous media is implemented for the theoretical 
analysis. The momentum balance theories proposed can explain 
pressure variation and pressure drop model can characterize the 
correlation between pressure variations and flow distribution. In 
addition, the influences of design parameters on flow distribu- 
tion regarding channel flow resistance, manifold width and air 
feed rate will be discussed in detail. 


2. Mathematical model 


Stack two-dimensional model is shown in Fig. 1(a) and (b). 
The gases enter manifolds through inlet ports, and then are dis- 
tributed into each cell and reach the electrodes through diffusion. 
After the electrochemical reactions are completed, un-reacted 
gases and products are discharged into manifolds from cells and 
then leave stacks through outlet ports. This flow field networks 
constitute an overall gas transport path. Besides simplifying 
stacks into two-dimensional model, other basic assumptions are 
as follows: 


the flow is steady; 

the flow is turbulent; 

the incompressible air is used as working fluid; 

electrochemistry, heat and mass transport phenomena are 

ignored; 

5. the gravity force is ignored and reference pressure for 1 atm 
is set; and 

6. the channels are filled with porous media. 


Pow PS 


Flow-field governing mechanism considered is turbulent. In 
order to simulate turbulent flow, k—e model is used to solve 
related transport equations. In addition, since channels are filled 
with porous media, treatment of porous media flow analysis is 
considered within turbulent model [20]. Governing equations 
are given as below: 
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Fig. 1. Schematic representation of fuel-cell stacks with manifold width (a) 
6mm, (b) 12mm. 


(b) Outlet | 


Continuity equation: 


V x (g=U) = 0 (1) 
Momentum equation: 
2 
U 
pV x UU) = —$¥p + Vx (on + SE" Q) 
ð oU; dU; 2 é6;=1 ifi=j 
T= fa( + i) - ody} : i 
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where « is turbulent kinetic energy; 
k (turbulent kinetic energy) equation: 


pV x (Uk) = V x (u it “+ o) + p- plde) 


2 k|U| 
+ ckp — (4) 

Oe 
where «k is turbulent kinetic energy, p= 
—pp(U; Ui, )(A(bU;)/ dx ;) is production rate of x, £ is turbulent 


dissipation rate, c,p¢*(k|U|//K) is additional production 
rate of « caused by porous media. o, = 1.0 and cy =0.28 are 
empirical constants. 

e (turbulent dissipation rate) equation: 


pV x (Us) = V x (u+ Hgo) $ Ca h 


EelpU) 2 
+ Cep (a — >) (5) 


where o, = 1.3, Cg, = 1.44, Ce2 = 1.92 are empirical constants. 

In the above governing equations, ¢ is porosity and K is per- 
meability of porous media. For the flow in the porous media, 
namely channel flow, ¢ is between 0 and 1 and K is a finite 
value depending on properties of porous media. For clear fluid 
flow, namely manifold flow, ġ is equal to one and K approaches 
infinitely large. At this stage, flow field will be entirely filled 
with fluid without the existence of porous media. 

The boundary conditions imposed to the boundaries of the 
computational domain are described as follows. The two air- 
feeding rates used are 100 and 300 standard liters per minute 
(SLPM). Air inlet velocity can be calculated from the feeding 
rate. Other boundary conditions are given as follows: 

Inlet boundary conditions: 

1. air average velocities at inlet ports are 3ms~! and 9ms7!, 
respectively; 
2. « and € can be defined as follows: 
3/2 

k= $ (UnTa?; e= a a=0.07L (6) 
where Uin is inlet velocity, Tin is turbulence intensity with 
value 0.1, L is hydraulic diameter with value 15 mm. 


Outlet boundary conditions are: 


. outlet pressure is constant at | atm; 

. fully developed flow at outlet ports, i.e. d6/dn=0, where n 
indicates normal direction relative to outlet ports, ® repre- 
sents physical properties such as velocities, pressure, etc. 


Ne 


Wall boundary conditions: 


1. flow is considered not to slip at wall therefore wall fluid 
velocity is zero; 

2. wall function is adopted to simulate turbulent wall shear 
stress. 


Since it is impossible to obtain an analytic solution of the com- 
plex convection-diffusion problem like this one, it will be solved 
by a finite volume method using a collocated cell-centered vari- 
able arrangement. The governing equations can be expressed in 
the form of a generalized transport equation 


V x (pug — IVb) = So (7) 


where @ is the general dependent variable, l'ẹ the exchange 
coefficient, Sg the source term, u velocity vector, and pis density. 
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With the discritization of the governing equations, the coupled 
finite-difference equations become 


app = Aude + aww + anon + asos + Sy (8) 


where @p is the value of ¢ at the current point P, dg ... ds stand 
for the values of the grid points adjacent to the point P, and ap 
... ag are known as the link coefficients. 

In this study, ESI-CFD ACE+ computational-fluid-dynamics 
software is used as modeling tools. The model is based on 
finite volume method, where the solution algorithm is SIMPLEC 
(Semi-Implicit Method for Pressure-Linked Equation Consis- 
tent) proposed by Van Doormal and Raithby [21] in 1984. 
Compared with SIMPLE (Semi-Implicit Method for Pressure- 
Linked Equation) developed by Patankar [22], SIMPLEC 
algorithm can provide quicker and more effective convergent 
solutions. 

Computational domains are divided into appropriate num- 
ber of meshes. Through this discretization process differential 
governing equations are transformed into algebraic equations 
on computational meshes. For computational accuracy and sta- 
bility, locally dense meshes should exist in specific domains 
with dramatically varied properties, such as channels, vari- 
able flow cross-sectional area or flow-field turns. According to 
grid tests, a body-fitted grid with 85000 computational cells is 
appropriate under compromise between modeling economy and 
reliability. 

In addition, continuity, momentum, turbulent kinetic and tur- 
bulent dissipation rate equations are solved in turn until the 
iterative process meets the convergence criteria. In this study, 
definition of convergence criteria indicates that the largest rel- 
ative error between consecutively two iterative residuals within 
overall computational domains is less than 1074. 


3. Results and discussions 


Using an advanced fuel-cell stack test station, overall output 
current, power, reactant flow rate, pressure drop, cell voltage, 
temperature distribution, etc. can be easily measured. As for 
reactant or coolant distribution, experimental device and tech- 
nique are generally difficult to implement so these microscopic 
phenomena have not yet been observed clearly. Some methods 
such as voltage distribution measurement can provide indirect 
information about reactant flow distribution. For example, while 
some cells in a stack produce lower voltage than others, it is pos- 
sible that less reactant is fed into the cells with low voltage. But 
other effects such as flooding or local overheating also might be 
the reason for the voltage loss. 

Physical mechanisms of flow distribution in stacks are dis- 
cussed and uniform degree of flow distribution is evaluated by 
changing design parameters. These parameters are included in 
channel flow resistance, manifold width and air feeding rate. 
Generally, channels with large flow resistance can cause large 
pressure drop and better performance in fuel cells. Channel flow 
resistance contributes mostly to pressure drop, but it is still 
not clear how channel flow resistance influences flow distribu- 
tion. In this study, different channel flow resistances caused by 


Table 1 

Stack model parameters and modeling conditions 

Variable Value 
Number of cells 72 

Gas channel depth (mm) 0.8 
Cell-to-cell distance in manifold axis (mm) 6 

Manifold width (mm) 6 and 12 
Width of inlet and outlet holes (mm) 15 

Porosity in cell channel 0.4 

Gas permeability in cell channel (m7) 2.5 x 10719 and 2.5 x 107° 
Back pressure (atm) 1 

Air flow rate (SLPM) 100 and 300 
Air viscosity (kg m7! s7!) 1.846 x 1075 
Air density (kg m~?) 1.1614 


adjusting permeability of porous media are analyzed to evalu- 
ate flow distribution. Manifolds are paths of reactant feed and 
discharge in fuel cell stacks. The geometrical size of manifold 
cross-sectional area will influence overall pressure drop in stacks 
but its effects on flow distribution is still not clearly known. In 
order to understand how flow distribution is affected by manifold 
width, different manifold sizes are used to evaluate this effect. In 
addition, this paper will also study the influence of different air 
feed on flow distribution. The model parameters and modeling 
conditions are given in Table |. The cases analyzed in this study 
are shown as below: 


Case 1: permeability K=2.5x107!°, manifold width 
D=6mm 

Case 2: permeability K=2.5x107!°, manifold width 
D=12 mm 

Case 3: permeability K =2.5 x 107°, manifold width D = 6 mm 
Case 4: permeability K=2.5 x 107?, manifold width 
D=12 mm 


Permeability of porous media can approximate the pressure drop 
in the channels. If the permeability is lower, pressure difference 
will be higher in the channels. Therefore, pressure drops caused 
by different channel configuration are realized when the per- 
meability is varied, i.e. the channel flow resistance is varied. 
By changing this parameter, the influence of channel flow resis- 
tance on flow distribution can be investigated. There is one order 
magnitude of difference for permeability in these cases so there 
will also be approximately one order magnitude of difference 
for the overall pressure drop. Larger deviation between magni- 
tudes of permeability can cause a more apparent difference of 
flow distribution. That is why we choose cases with permeabil- 
ity K=2.5 x 107!° and K=2.5 x 107°. In addition, this study 
also evaluates how different manifold widths influence flow dis- 
tribution. Manifold widths of 6mm and 12 mm are evaluated 
and 6mm is the origin design in stacks. In order to observe 
the apparent effect of manifold width to flow distribution, cases 
where size of origin width are doubled, i.e. 12 mm, are investi- 
gated. In the following section, flow pressure variation, channel 
pressure drop and cell flow rate under 300 SLPM air feeding 
rate will be discussed separately. 
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3.1. Flow pressure in stacks 


Flow pressure variations of case 1, case 2, case 3 and case 
4 under 300 SLPM air feeding rate are shown in Figs. 2—4 and 
Fig. 5(a)-(c), respectively where (a) is overall flow pressure in 
stacks and (b) and (c) are respectively average pressure variation 
along inlet and outlet manifold length. 

For overall flow pressure in stacks, i.e. the right hand side 
of Figs. 2-5, the highest pressure occurs in the inlet manifold 
and the lowest one in the outlet manifold. In channels, pres- 
sure gradually decreases from the inlet manifold to the outlet 
one. On comparing the overall pressure drop of each case, case 
1>case 2>case 3>case 4. This shows that the increased chan- 
nel flow resistance due to lower permeability contributes to a 
larger overall pressure drop. As mentioned earlier, the pres- 
sure drop in stacks is mainly controlled by the permeability i.e. 
K=2.5 x 10~!° (case 1 and case 2) and K=2.5 x 107° (case 3 
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and case 4). These results show that one order decrease of mag- 
nitude for permeability will increase the overall pressure drop 
by about one order magnitude, i.e. pressure drop increases from 
3377 Pa (case 4) and 4003 Pa (case 3) to 29540 Pa (case 2) and 
30160 Pa (case 1). 

Larger manifold widths cause less overall pressure drop. 
While manifold width increases from 6mm (case 1 and case 
2) to 12 mm (case 3 and case 4), the overall pressure decreases 
from 4003 Pa (case 3) and 30160 Pa (case 1) to 3377 Pa (case 
4) and 29540 Pa (case 2). According to fluid mechanics theo- 
ries, lowered fluid velocity due to enhanced manifold width will 
contribute to less pressure drop. 

On the left hand side of the figures, the top plot is pres- 
sure variation along the inlet manifold and the bottom one is 
pressure variation along the outlet manifold. The x-coordinate 
indicates distance along manifold length where the origin is 
at the feeding end of stacks, and the y-coordinate is pressure 


P — N/m*2 
3.016E+004 


3£+00 


2.5E+004 


2E+004 


1.5£+004 


1E+004 


5000 


-73.75 


Oi 


(Pa) f 


a 250 
Pressure variations along 
outlet manifold length 


s0 a it L 4 i L L 4 
0 005 Ol 015 02 025 03 035 04 045 


Chordlength (m) : 


Fig. 2. (a) Stack flow pressure, (b) pressure variations along inlet manifold length, (c) pressure variations along outlet manifold length for case 1 (K=2.5 x 107!°, 


D=6mm, 300 SLPM). 
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Fig. 3. (a) Stack flow pressure, (b) pressure variations along inlet manifold length, (c) pressure variations along outlet manifold length for case 2 (K=2.5 x 107!°, 


D=12mm, 300 SLPM). 


with unit: Pa (N m°). Pressure increases from the feeding end 
to the closed one of stacks at either inlet or outlet manifolds. 
Gas flow rate gradually decreases along inlet manifold length 
due to gas feeding into each cell. The decreased flow rate, i.e. 
decreased fluid velocity will cause a pressure gradient opposite 
to flow direction. For the outlet manifold, gas flow rate gradually 
increases from the closed end to the discharged one due to gas 
discharge from each cell. The increased flow rate, i.e. increased 
fluid velocity will cause a pressure gradient parallel to flow 
direction. 

The pressure difference in manifolds is mainly associated 
with manifold width and is less dependent on the permeability, 
i.e. channel flow resistance. Larger manifold width causes a less 
pressure difference along the manifold length. By comparing 
cases with a width of 6mm (case | and case 3) to cases with 
a width of 12mm (case 2 and case 4), it shows that the pres- 
sure difference along the inlet manifold decreases from 140 Pa 
(case | and case 3) to 40 Pa (case 2 and case 4), and that along 
the outlet manifold from 550 Pa (case | and case 3) to 130 Pa 


(case 2 and case 4). This is because wider manifolds cause 
lower fluid velocity which leads to less pressure difference in 
manifolds. 

The above analysis of manifold pressure difference shows 
that the outlet manifold contributes to a larger pressure difference 
than in the inlet one. This will cause an uneven pressure-drop in 
the cells, which will lead to uneven gas fed in each cell, i.e. an 
uneven flow distribution. 


3.2. Momentum balance theories 


Momentum balance theories [23] characterize overall 
momentum balance relationship consisting of momentum flux, 
pressure, and wall friction in a specific control volume. In the 
theories, net momentum flow rate must be equal to the sums of 
all forces, and the governing equation is: 


Div 
1 PAV y= YÈ (9) 
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Fig. 4. (a) Stack flow pressure, (b) pressure variations along inlet manifold length, (c) pressure variations along outlet manifold length for case 3 (K=2.5 x 107°, 


D=6mm, 300 SLPM). 


where m is mass flow rate, V is fluid velocity, and 5> F is the 
net value of all forces. 

For steady flows, a momentum balance exists between net 
momentum flow rate, net forces of pressure and wall shear stress 
as follows: 


SomvV=S)°-PS+ YO 


where P is pressure, S is pressure action area, and Fe is wall 
shear stress. 

Inlet and outlet manifolds are chosen as control volumes, 
as shown in Figs. 6 and 7, respectively. The two momentum- 
balance relationships are: 


(10) 


—mV = (Pi — Po)S — Fy = —AP,S — Fry 
mV — Fry 


> AP, = — 


a1) 


mV = (Py — P3)S — Fp = AP S — Fp 
mV + Fp 


> AP = 5 


(12) 


In the two equations, rn V represents feeding momentum flow 
rate at the inlet manifold or discharging momentum flow rate 
at the outlet one. The Fs; and Fp are friction force generated 
by fluid flow through inlet and outlet manifolds. The Pı and 
P» are pressures at two ends of the inlet manifold, and P3 and 
P4 are pressures at two ends of the outlet one. The S is pres- 
sure action area. According to the momentum balance theories, 
the deviation between pressure difference in the inlet manifold: 
AP, =P 2 — P, and that in the outlet manifold: AP, = P4 — P3 
is: 


Fry + Fr 


AP, = AB, — < AP, (13) 
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Fig. 5. (a) Stack flow pressure, (b) pressure variations along inlet manifold length, (c) pressure variations along outlet manifold length for case 4 (K=2.5 x 107°, 


D=12mm, 300 SLPM). 


This theory could explain the higher pressure drop of the 
outlet manifold than the inlet manifold. Causes for the uneven 
pressure drop among the cells also could be explained. 


3.3. Cell pressure drop and mass flow rate 


In gas-distributed processes, forces caused by pressure gradi- 
ents in channels drive the gas into each cell from the manifolds. 
In order to confirm that cell pressure drop dominates flow dis- 
tribution, comparison between cell pressure drop and mass flow 
rate was studied. 


Closed end 


Fig. 6. Schematic representation of inlet manifold control volume. 


Fig. 8 is a plot of cell pressure drop under 300 SLPM air 
feeding rate. A dimensionless pressure drop is implemented to 
progress comparisons between each case. It can be defined as 
follows: 
AP AP AP (14) 

~ AP SS AP/N 
where AP is dimensionless pressure drop, AP is cell pressure 
drop, AP = X` AP/N is average cell pressure drop, and N is 
number of cells. 


Closed end Outlet 


Fig. 7. Schematic representation of outlet manifold control volume. 
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Fig. 8. Dimensionless cell pressure drop for each case (300 SLPM). 


In plots of cell pressure drop, the x-coordinate indicates the 
cell number, where no. 1 is indicated as the cell at the feed- 
ing end and no. 72 as the cell at the closed one in stacks. 
The y-coordinate is defined as dimensionless pressure drop. The 
magnitude of dimensionless pressure drop for different cases are 
case 1: 1.009-0.995, case 2: 1.002—0.999, case 3: 1.081—0.960, 
case 4: 1.018—0.991. Dimensionless pressure-drop variation of 
all cases gradually decreases from the feeding end to the closed 
one in stacks. As for the uniform degree of pressure-drop distri- 
bution, the best one is case 2, and the next ones are case | and 
case 4, and the worst one is case 3. 

The dimensionless mass flow rate in the individual cell for the 
300 SLPM air feeding rate is given in Fig. 9. A dimensionless 
mass flow rate is implemented to compare the flow through the 
individual cells. It can be defined as follows: 

e m m 

is s (15) 
m Q/N 

where M is dimensionless mass flow rate, rn is cell mass flow 

rate, and rn = Q/N is average cell mass flow rate, where Q is 

total mass flow rate of stacks and N is number of cells. 
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Fig. 9. Dimensionless cell mass flow rate for each case (300 SLPM). 


In the plot of cell mass flow rate, the x-coordinate indicates 
the cell number, where no. 1 is the cell at the feeding end and 
no. 72 as the cell at the closed one in the stack. The y-coordinate 
is defined as dimensionless mass flow rate. The magnitude of 
dimensionless mass flow rate for the different cases are: case 1: 
1.008—0.996, case 2: 1.002—0.999, case 3: 1.068—0.967 and case 
4: 1.015-0.992. Variation of dimensionless mass flow rate for all 
cases gradually decreases from the feeding end to the closed one 
in stacks. Uniform degree of flow distribution is almost the same 
as that of pressure-drop variation. This confirms that the driving 
forces governing gas feed into each cell results from pressure 
drop in channels. 


3.4. Analysis of flow distribution 


A pressure drop model shown in Fig. 10 is constructed to 
investigate pressure variation where manifolds, channels and gas 
flow direction are indicated, respectively. Since flow pressure is 
defined relative to the outlet pressure, pressure at outlet ports 
is set to zero. (The red solid line indicates pressure drop of 
cell no. 1, i.e. AP}; the red dotted line indicates pressure drop 
of cell no. n, i.e. AP,, the green solid line indicates pressure 
difference along the inlet manifold length, i.e. AP,; and the blue 
solid line indicates pressure difference along the outlet manifold 
length, i.e. AP,.) The following relationship will be satisfied for 
pressure drop AP, AP;, AP, and APh: 


AP, = AP; + AP, — APh (16) 


The above equation can evaluate qualitatively the effect of chan- 
nel flow resistance and manifold width on flow distribution. 
The ratio of AP, to AP; is calculated by using the relation- 
ship: AP, = AP; + AP, — AP» from pressure-drop model and 
AP, < AP» from momentum balance theory respectively: 

AP, AP, +AP,— Ah AP, — AP, 


= =14 22b 1 17 
AP, AP 7 AB a) 


The above equation defines the ratio of pressure drop of cell 
no. n to no. 1. This ratio is less than one, i.e. cell pressure drop 
decreases gradually from the feeding end to the closed one, and 
it also corresponds with the modeling results. 

From the above discussion, it can be seen that the cell pres- 
sure drop dominates gas flow into each cell. Therefore, the ratio 


manifold qm Gas flow directio = 


Fig. 10. Schematic representation of stack pressure drop model. 
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Table 2 
Stack design parameters and corresponding pressure drop (300 SLPM) 
Permeability Manifold width 
D=6mm D=12 mm 
AP, (Pa) 139 AP, (Pa) 38 
K=2.5 x 107! m? AP, (Pa) 550 AP, (Pa) 130 
AP; (Pa) 29691 AP; (Pa) 29476 
AP, (Pa) 140 AP, (Pa) 39 
K=2.5x 107° m? AP, (Pa) 535 AP, (Pa) 128 


AP; (Pa) 3530 AP; (Pa) 3316 


between cell pressure drops can be qualitatively used to evaluate 
gas flow distribution. The ratio of mass flow rate of cell no. n to 
no. 1 can be defined as follows: 


itn — APn AP, — AP, 
mo AP, AP; 


(18) 


The above equation indicates the ratio of cell mass flow rate 
between the two cells which can be calculated with manifold 
pressure differences AP, and AP» and cell pressure drop AP}. 
When this ratio is closer to one a more uniform flow distribu- 
tion will occur in stacks. Next, this equation is implemented to 
analyze the above cases. 

Model design parameters and contributed pressure drops 
under 300 SLPM air feeding rate are shown in Table 2. Pressure 
differences in the inlet and outlet manifold, i.e. AP, and APy 
are dominated by manifold widths. Compared with manifold 
widths, permeability, i.e. channel flow resistance has a weaker 
influence on manifold pressure difference but dominates cell 
pressure drop. A decrease of one order magnitude in perme- 
ability will result an increase of one order magnitude of overall 
pressure drop. But variable manifold widths contribute almost 
nothing to overall pressure drop. The ratios of cell mass flow 
rate for each case under 300 SLPM air feeding rate are shown 
in Table 3 and no. n is defined as 72, where case | is 0.986, 
case 2 is 0.997, case 3 is 0.888, and case 4 is 0.973. The ratios 
of each case correspond qualitatively with the modeling results. 
As for the uniform degree of flow distribution, case 2 is the best 
one followed by case 1 and case 4 in turn and case 3 is the 
worst one. 

From the above analysis, low permeability can contribute 
to a large cell pressure drop AP; and improve flow distribution. 
Large manifold width can cause less deviation between inlet and 
outlet manifold pressure difference, i.e. APp — AP, and then it 
also will improve flow distribution. 


3.5. Flow distribution under low air feeding rate 


The above cases are analyzed under 300 SLPM air feeding 
rate and the stack will generate 4-5 kW output power at these 
feeding rate. In order to study stack flow distribution for output 
power 1-2 kW, approximately 100 SLPM air feeding rate was 
selected for investigation. The inlet air velocity was 3ms7! at 
100 SLPM air feeding rate. Besides the air feeding rate, stack 
model parameters and modeling conditions are the same as the 
above cases. In this section, the difference of flow distribution 
under different air feeding is studied by momentum balance the- 
ories and the influences of channel resistance and manifold width 
on flow distribution are analyzed. 

Flow pressure variation for case 1, case 2, case 3 and 
case 4 under 100 SLPM air feeding rate are shown in 
Figs. 11—-14(a)-(c), respectively where (a) is overall flow pres- 
sure in stacks and (b) and (c) are respectively average pressure 
variation along inlet and outlet manifold length. 

Under 100 SLPM air feeding, overall pressure drop of dif- 
ferent cases is case 1 >case 2>case 3>case 4. This is because 
larger channel flow resistance contributes to a larger overall pres- 
sure drop. While manifold width is enhanced lower fluid velocity 
contributes to less pressure drop. 

Pressure increases from the feeding end to the closed one in 
stacks at either inlet or outlet manifolds. Decreased fluid velocity 
from the feeding end to the closed one along the inlet manifold 
length will cause a pressure gradient opposite to flow direction. 
For the outlet manifold, increased fluid velocity from the closed 
end to the feeding one will cause a pressure gradient parallel to 
velocity direction. 

Overall pressure drops of each case under 100 and 300 SLPM 
air feeding rate are shown in Table 4. 

Since channel flow resistance, i.e. permeability contributes 
largely to flow pressure drop, the correlation between pressure 
gradient, flow velocity and permeability can be approximated 
with the equation: Vp% @uU/K. From this equation, overall 
pressure drop is roughly proportional to air feed, i.e. the ratio 
of overall pressure drop for 300 SLPM to that for 100 SLPM is 
about 3. For realistic channels, such as parallel, serpentine, or 
Z pattern channel pressure drop is also roughly proportional to 
flow velocity. ‘, channel flow resistance approximated by per- 
meability is satisfied qualitatively with governing mechanism of 
channel pressure drop. 

Fig. 15 is a plot of cell pressure drop under 100 SLPM air 
feeding rate. A dimensionless pressure drop is implemented to 
progress comparisons between each case. Cell pressure drops of 
each case under 100 and 300 SLPM air feeding rate are shown 
in Table 5. 


Table 3 Table 4 

The ratio of cell mass flow rate for each case (300 SLPM) Overall pressure drops of each case under 100 and 300 SLPM air feeding rate 
Case Ratio Air feeding rate 100 SLPM (Pa) 300 SLPM (Pa) 
1 0.986 Case 1 9468 30160 

2 0.997 Case 2 9430 29540 

3 0.888 Case 3 1090 4003 

4 0.973 Case 4 1052 3377 
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Fig. 11. (a) Stack flow pressure, (b) pressure variations along inlet manifold length, (c) pressure variations along outlet manifold length for case 1 (K=2.5 x 107!°, 


D=6 mm, 100 SLPM). 


Under different air feed dimensionless pressure-drop varia- 
tion of all cases gradually decrease from the feeding end to the 
closed one in stacks. As for the uniform degree of pressure-drop 
distribution, the best one is case 2 and the next are case 1 and 
case 4 in turn, and the worst one is case 3. In addition, lower air 
feed will cause more uniform cell pressure-drop variation, and 
this will also cause more uniform flow distribution. 

Fig. 16 is a plot of cell mass flow rate under 100 SLPM air 
feeding rate. A dimensionless cell mass flow rate is implemented 


to study comparisons between each case. Cell mass flow rate of 
each case under 100 and 300 SLPM air feeding rate are shown 
in Table 6. 

Under different air feed variation of dimensionless mass flow 
rate of all cases gradually decreases from the feeding end to 
the closed one in stacks. Uniform degree of flow distribution 
is almost the same as that of pressure-drop distribution and 
lower air feed indeed causes more uniform distribution of cell 
flux. 


Table 5 Table 6 

Cell pressure drops of each case under 100 and 300S LPM air feeding rate Cell mass flow rate of each case under 100 and 300 SLPM air feeding rate 
Air feeding rate 100 SLPM 300 SLPM Air feeding rate 100 SLPM 300 SLPM 
Case 1 1.004—0.998 1.009-0.995 Case 1 1.004—0.998 1.008-0.996 
Case 2 1.0008-0.9996 1.002-0.999 Case 2 1.0008-0.9996 1.002-0.999 
Case 3 1.036-0.983 1.08 1-0.960 Case 3 1.034-0.984 1.068-0.967 
Case 4 1.007-0.996 1.018-0.991 Case 4 1.007-0.997 1.015-0.992 
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Fig. 12. (a) Stack flow pressure, (b) pressure variations along inlet manifold length, (c) pressure variations along outlet manifold length for case 2 (K=2.5 x 107!°, 


D=12mm, 100 SLPM). 


Momentum balance theories are used to study the difference 
of flow distribution under different air feed and the influence of 
channel resistance and manifold width on flow distribution. The 
ratio of mass flow rate of cell no. n to no. 1 can be defined with 
manifold pressure differences AP, and AP, and cell pressure 
drop AP as follows: 


AP, — AP, 
z1 (19) 
AP 


Mn ns AP, 
mı AP, 


Channel pressure drop AP is approximately proportional to 
flow velocity, i.e. AP; «x V. However, manifold pressure differ- 
ences AP, and AP» are roughly proportional to the square of 
flow velocity, i.e. AP, or APp « V2. This is because manifold 
pressure differences AP, = mV — Fy, /S and AP = mV — 
F ¢2/S can be derived respectively as follows: 


mV — Fa _ pV? -Ek V’/2 


AP, = x y? 20 
a 5 5 x (20) 


mV + Fe pV +Y kV?/2 


y? 21 
5 5 (21) 


AP, = 


where momentum flux is proportional to the square of flow 
velocity. According to fluid mechanics theories [24], friction 
forces are approximated as multiplication of friction coefficient 
and the square of flow velocity, i.e. Fs kV7/2, where k is fric- 
tion coefficient. Therefore, the manifold pressure difference is 
also proportional to the square of flow velocity. 

Model design parameters and contributed pressure drop 
under 100 SLPM air feeding rate are shown in Table 7. By com- 
paring Table 4 to Table 2, cell pressure drop under 300 SLPM air 
feeding rate is three times larger than that under 100 SLPM one, 
and manifold pressure difference under 300 SLPM air feeding 
rate is nine times larger than that under 100 SLPM one. This 
corresponds to conclusions derived from momentum balance 
theories. While the air feed is varied and the other parameters are 


+. in y APn ~ q _ APAP, ; 
fixed, the equation: m REC 1 AP; can be described 
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Fig. 13. (a) Stack flow pressure, (b) pressure variations along inlet manifold length, (c) pressure variations along outlet manifold length for case 3 (K=2.5 x 107°, 


D=6 mm, 100 SLPM). 


with a flow velocity function as follows: 


fitn 2  APo— APs ACV) IV) 
mı AP, fV) ' 
where f(V)& V (22) 
Table 7 
Stack design parameters and corresponding pressure drop (100 SLPM) 
Permeability Manifold width 
D=6mm D=12 mm 
AP, (Pa) 12 AP, (Pa) 4 
K=2.5 x 107! m? AP» (Pa) 69 AP, (Pa) 16 
AP; (Pa) 9456 AP; (Pa) 9425 
AP, (Pa) 12 AP, (Pa) 4 
K=2.5 x 107! m? AP, (Pa) 68 AP, (Pa) 15 
AP; (Pa) 1078 AP, (Pa) 1048 


When the air feed is reduced, decreased flow velocity makes 
J(V) smaller to lead the ratio: 71,,/m, to approach close to one. 
By comparing Table 5 to Table 3, flow distribution under 100 
SLPM air-feeding rate is more uniform. The ratios of cell mass 
flow rate for each case in Table 8 are, case 1: 0.986, case 2: 
0.997, case 3: 0.888, and case 4: 0.973. The ratios qualitatively 
correspond to variations of cell mass flow rate in Fig. 16. As for 
each design parameter, low permeability contributes to a large 
cell pressure drop AP, and improves flow distribution. Large 


Table 8 

The ratio of cell mass flow rate for each case (100 SLPM) 
Case Ratio 
1 0.994 
2 0.999 
3 0.948 
4 0.990 
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Fig. 14. (a) Stack flow pressure, (b) pressure variations along inlet manifold length, (c) pressure variations along outlet manifold length for case 4 (K=2.5 x 107°, 
D=12mm, 100 SLPM). 
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Fig. 15. Dimensionless cell pressure drop for each case (100 SLPM). Fig. 16. Dimensionless cell mass flow rate for each case (100 SLPM). 


C.-H. Chen et al. / Journal of Power Sources 173 (2007) 249-263 263 


manifold width causes a less deviation between the inlet and 
outlet manifold pressure difference, i.e. AP) — AP, and then 
also improves flow distribution. 


4. Conclusions 


A fuel-cell stack model is constructed using computational 
fluid dynamics to investigate pressure variation, flow distribu- 
tion for different channel flow resistance, manifold widths and 
air feed. Proposed momentum balance theories and pressure- 
drop model can be used to qualitatively investigate the flow 
distribution occurred in stacks. 

Following are the conclusions of this study.This study eval- 
uates flow distribution with a pressure-drop model and explains 
the physical mechanism governing this phenomenon qualita- 
tively.Different air feed will cause different flow distribution. 
This study indicates that lesser air feed promotes more uniform 
flow distribution than higher air feed. This is because there exist 
different governing mechanisms between manifold pressure dif- 
ference and channel pressure drop from momentum balance 
theories and pressure-drop model.Different channel design in 
stacks will affect flow distribution. Channels with large flow 
resistance contribute more pressure drop, and then cause a more 
uniform flow distribution. Therefore, channel design with large 
flow resistance is advantageous for flow distribution. While man- 
ifold widths increase, a more uniform flow distribution will be 
achieved. For manifold design in stacks, its width should be 
enhanced as much as possible and this is also beneficial for 
lowering overall stack pressure drop.The important points to be 
considered in stack design are lowering overall pressure drop and 
attaining high performance and stable operation. However, com- 
petitions and conflict always exist between these performance 


characteristics. In order to satisfy multiphase requirements of 
stack design, enhanced manifold size is a better solution to pro- 
mote flow distribution and implement flow-field optimal design 
of fuel-cell stacks. 
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